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Proponents of Complexity Science believe that the huge variety of emergent phenomena observed 
throughout nature, are generated by relatively few microscopic mechanisms. Skeptics however point 
to the lack of concrete examples in which a single mechanistic model manages to capture relevant 
macroscopic and microscopic properties for two or more distinct systems operating across radically 
different length and time scales. Here we show how a single complexity model built around cluster 
coalescence and fragmentation, can cross the fundamental divide between many-body quantum 
physics and social science. It simultaneously (i) explains a mysterious recent finding of Fratini et al. 
concerning quantum many-body effects in cuprate superconductors (i.e. scale of 10 -9 — 10 -4 meters 
and 10 -12 — 10 -6 seconds), (ii) explains the apparent universality of the casualty distributions in 
distinct human insurgencies and terrorism (i.e. scale of 10 3 — 10 6 meters and 10 4 — 10 8 seconds), 
(iii) shows consistency with various established empirical facts for financial markets, neurons and 
human gangs and (iv) makes microscopic sense for each application. Our findings also suggest 
that a potentially productive shift can be made in Complexity research toward the identification of 
equivalent many-body dynamics in both classical and quantum regimes. 

The central limit theorem predicts that a collection of objects in which there are no hidden correlations (e.g. 
a handful of coins being tossed) will produce fluctuations which have an approximate Gaussian distribution. By 
contrast, the fluctuations emerging from systems containing correlations which cross multiple length and/or time 
scales can exhibit significant deviations from Gaussian behavior. A statistical form which emerges from many distinct 
systems is the exponentially truncated power-law distribution p{x) oc x~ a ex.p(— x/xq) where a typically takes values 
between one and foui?™. It is tempting to infer that if two such systems A and B share similar values a a and 
then they also share a common underlying dynamical mechanism. But this is not generally true. Indeed one of the 
most damaging criticisms of Complexity Science as a unifying discipline, is that mechanisms which make sense in the 
context of system A may make little sense for system B. For example, a two-dimensional model with nearest-neighbor 
interactions seems reasonable for vehicle traffic or fish shoals, but not for herding in global financial markets where 
interactions become essentially independent of spatial separation^L^. 

Complex phenomena appear across all length and time scales. In the physics community, the toughest complexity 
arguably lies at the level of quantum many-body phenomena, such as high temperature superconductivity in the 
cuprates^l i n the social sciences, it arguably lies in the field of human conflict 6 . A concrete quantitative connection 
between these two phenomena would seem highly unlikely - yet it is precisely this connection that we uncover here. 
Specifically, we use a single dynamical model to simultaneously provide a first explanation for the remarkable yet 
mysterious recent finding of Fratini et a/.Hl (see Fig. 1(a)) concerning the complex, nanoscale quantum- mechanical 
world of high-temperature superconductivity, and an explanation of a recent empirical observation concerning casualty 
distributions in human insurgencies and terrorism 3 . In addition to making microscopic sense for each system in 
turn, the model successfully reproduces the distributions reported for several other complex phenomena of current 
interest^™ 

X-ray diffraction (XRD) studies by Fratini et al. 11 in the high-temperature superconducting (SC) cuprate system 
La2Cu04+ 2y revealed a partial ordering of a stripe-like L Q2 superstructure' of the interstitial oxygen atoms (i-Os), 
depending on the heat treatment of the samples. The probability distribution of the Q2 XRD intensity, x = I(Q2)/Iq 
(normalized to the intensity Io of the tail of the main crystalline reflections at each spatial point), was shown 11 
to exhibit fractal scaling over a wide range: p(x) oc x _a exp(— x/xq) with a = 2.6 ± 0.2 (see Fig. 1(a)). Since 
the intensity is proportional to the volume of scatterers, this suggests that the distribution of nanoscale volumes of 
ordered i-Os follows the same power-law pattern. Thus x can be re-interpreted as the volume occupied by Q2-ordered 
i-Os. Remarkably, it was found 11 that this fractal arrangement of Q2-ordered i-Os promotes superconductivity in 
La2Cu04 +2y (i.e. the SC transition temperature T c increases when the fractality is more prominent). 

Figure 1(b) shows the prediction of the complexity model which is solved analytically in Appendix A, while the 
common underlying mechanism is shown schematically in Fig. 2. The field of high-T c superconductivity is packed with 
a variety of theories^, ranging from the conventional BCS approach through to the suggestion of behavior associated 
with special black holes 13 . However, as stated in the accompanying description 13 to Fratini et aZ.'s paper^, 'there is 
nothing in the textbooks even hinting at an explanation' for why fractal-defect structure, ranging from a micrometer 
up to fractions of a millimeter scale, should cause the observed enhancement of T c - and in particular, why the valued 
a = 2.6 ± 0.2 (Fig. 1(a)) emerges. 
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FIG. 1: (a) Experimental results (triangles and squares) compared to an exponentially truncated power-law curve (black line) 
having a = 2.6 db 0.2. Triangles correspond to T c = 40 K and squares correspond to T c = 16 + 32 K. Information extracted 
from Fig. 2c of Fratini et a/.= and shows data points lying above the lower cutoff, (b) Analytic solution of the model (Eq. 
(1)) corresponds to an exponentially truncated power- law with a — 2.5 (see Appendix A for full discussion). For illustration 
purposes, (b) shows z/f rag = v = 1 — zy CO ai but the same form emerges for more general values. Following Fratini et aZ.™ the 
data in (a) are scaled by xo and the full exponential tail is not shown. By contrast, (b) shows the full unsealed theoretical 
form. 



Experimental evidence^"^ indicates that cuprate superconductors are characterized by a temperature range {Tfl uc }, 
above T c , where SC fluctuations exist but the global SC phase coherence is incomplete, and that the SC-fluctuation 
regime does not coincide with the pseudogap phase of the cuprates which probably reflects a two-gap scenario^ 
(see discussion elsewhere^). Furthermore, it has been demonstrated that effects such as compressive strain on thin 
films^l can raise T c into the SC-fluctuation regime. This regime is therefore characterized by the existence of Cooper 
pairs (CPs) which lack complete phase coherence throughout the system, with the number of CPs depending on 
temperature (T) in the thermodynamic limit. The ordering of the i-Os in La2Cu04 +? J^has a stabilizing effect on the 
coherence between the CPs which is diminished by the existence of disorder; furthermore, the location of the ordered 
i-Os between the Cu02 planes is helpful for the establishment of inter-planar coherence. Hence an individual ordered 
mesoscopic i-0 region will provide a favorable environment for the portion of the SC state which finds itself within it, 
implying that the sample-wide many body state will have spatially heterogeneous coherence depending on how the 
ordered i-0 regions are arranged. Our challenge is therefore to explain how a fractal space with a = 2.6 ±0.2 manages 
to protect the many-body phase coherence (and therefore long-range correlations) at temperatures for which an SC 
state would ordinarily be very fragile. 

Irrespective of the actual microscopic pairing mechanism giving rise to coherent CPs, one can therefore model the 
{^fluc} regime as a dynamical competition between (1) the entropy- driven tendency of a coherent SC state to break 
up (i.e. fragment) through phase-breaking events, into smaller clusters where the intra-cluster phase coherence is 
strong but the inter-cluster coherence is weak; and (2) the tendency of the underlying SC mechanism to synchronize 
the phases within two clusters, forming one larger coherent droplet as shown schematically in Fig. 2. Close to T c , 
(2) will dominate, while considerably above it (1) will dominate - however both tendencies coexist for T £ {Tfl uc }. 
Hence a T-dependent population of N ^> 1 CPs comprises, for T G {Tfl uc }, a set of clusters with each one having 
its own internal phase coherence. Each cluster has a size s representing the number of phase-coherent objects (i.e. 
CPs), and the number n s (t) of such clusters evolves over time as a result of fragmentation and coalescence (i.e. (1) 
and (2) above) such that ^Zf = i sn s (t) = N at each timestep. At each timestep, a coherence-breaking event will 
occur with probability ^f rag in a cluster which is randomly chosen according to its size s, mimicking the fact that 
larger clusters will suffer proportionally more quantum fluctuations. In short, a cluster with more members has more 
chances of initiating a dephasing event. With probability ^ CO ab this cluster instead coalesces with another one chosen 
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FIG. 2: Schematic of model processes of coalescence and fragmentation by which coherence of two clusters becomes synchronized 
(i.e. coalescence) or the coherence within a single cluster is suddenly destroyed by a decohering event (i.e. fragmentation). 
These processes are generic, but the language chosen is for specific application to superconductivity. 



randomly according to its size. This mimics the tendency of two clusters to synchronize their phase coherence as the 
system heads toward global phase-coherence in the ideal SC state, i.e. any subsequent coalescence events will likely be 
initiated by pairwise coherence between individual members in the two groups and hence the probability will depend 
on the number of members. We stress that we use the term 'coalescence' to simply mean that two groups act in a 
coordinated way, not necessarily that they are physically joined. As T is lowered towards T c , ^f rag is likely to decrease 
due to the reduced frequency of coherence-breaking processes, while v coa \ is likely to increase as the establishment of 
a coherent SC state is approached. 

The Master Equations for this model are: 

f) s— 1 „ oo 

OU^ _ Z^oal ST^u ( k) Vfr a gSn s ZV coal STl s ^ 

Q t ~ TV 2 Z^ KUk ^ S K ) n s-k N N 2 Z^ Knk ' S - A ' W 

k=l k=l 

r\ oo oo 

~dt ~ 1^2^ knk -^^2^ knk > 

k=2 k=l 

where we make the sensible approximation that TV — >> oo. Simpler versions of Eq. (1) have been presented before^^, 
however this is the first treatment with general values for ^f rag and v C03 \ and it is the first application to superconduc- 
tivity. The distance-independence is justified by the fact that the phase coherence in a macroscopic SC state can exist 
over large distances (e.g. 10 _9 m-10 _4 m) even if it is fragile (Fig. 2), and hence makes Eq. (1) equivalent to a mean- 
field model (i.e. no distance dependence). Solving the model analytically in the steady-state limit (i.e. dn s /dt = 0) 
yields a closed-form expression for the time-averaged probability distribution for phase-coherent cluster sizes which 
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FIG. 3: Model prediction (Eq. (2)) for how the exponential cutoff so, the fragmentation probability z/f rag and the coalescence 
probability ^coai are related. A portion of the possible parameter value space is shown for illustration. 

is in exact agreement with Fratini et a/. 11 , i.e. it yields an exponentially truncated power-law p(s) oc s~ a ex.p(— s/sq) 
with a = 5/2 = 2.5 (see Appendix A). We also uncover the following novel theoretical relationship between cutoff sq 
and fragmentation and coalescence probabilities ^f rag and v coa \: 



so = ~ 



In 



4(Vfrag + ^coal)^coal 
(^frag + 2^ coa i) 2 



(2) 



The model's 5/2 truncated power-law form is remarkably robust to generalizations, e.g. coalescence of multiple 
groups, fragmentation into groups larger than one, a slowly time-varying particle number iV, and it holds for a wide 
range of ^f rag and ^ coa i value d 2 W 24 * . It does not depend on the size of the individual fragmented parts, as long as they 
are all small, since the size of the largest pieces just dictates the value of s above which the 5/2 power-law result kicks in. 
Figure 1(b) illustrates the model's predicted analytic form, for the simple demonstrative case of ^f rag = v = 1 — ^ C oai- 
As ^f rag increases and/or v C03 \ decreases, the cutoff sq sets in at lower s and eventually dominates the power-law. (For 
Fig. 1(a), adapted from Fratini et a/. 11 , the axis rescaling shifts the red experimental curve to the right and the tails 
fall outside the range shown, as opposed to Fig. 1(b) which shows the full unsealed theoretical form including the 
exponential tail). Equation (2) is illustrated in Fig. 3. In addition to reproducing Fratini et al 'J^ precise functional 
form, the model allows us to interpret the cutoff sq in terms of the dynamic processes of fragmentation and coalescence 
(see Fig. 3). 
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Our theory does not require any specific behavior at the level of individual CPs, nor precise details or mechanisms 
involving the location of the i-Os between the Cu02 planes. Instead, it takes the novel approach of simultaneously 
describing the partially coherent SC state across several orders of lengthscale magnitude, predicting that it can fit 
perfectly into the fractal space created by the ordered i-Os such that each cluster resides within a pocket full of 
ordered i-Os, as sketched in Fig. 2. Since the fractal space of the ordered i-Os (which develops in the samples at high 
temperatures where the heat treatment is applied^) is basically static for T G {Tfl uc }, the coupling of the fractal 
space of the CPs to it results in an increase of SC coherence (see above). Consequently T c is raised within the {Tfl uc } 
temperature regime, as has been observed^. In samples where the fractal i-0 box is very incomplete, T c is not raised 
throughout an entire sample which hence results in the observed mixed stated 

A question to be clarified in future research is whether the ordered i-0 cluster formation process follows a similar 
mechanism to that of the CP cluster formation. Similarly to the CPs, the collective (phonon-like) excitations, 
associated with the ordering of the i-Os, have a quantum nature and can in principle exhibit long-range phase 
coherence. The mobility of the i-Os is quite high at the heat-treatment temperatures^, and fluctuations of the Q2- 
ordering phase transition occur. However, the ordering phase transition cannot be completed throughout the sample 
since the mobility is reduced when T is lowered - hence an i-0 fractal structure gets locked in at low T. 

The Master Equations in Eq. (1) also provide a plausible explanation for the recently observed tendency of modern 
insurgent groups to produce approximate power-law casualty distributions witfP a « 2.5. In this human context, 
the distance-independent coalescence of Eq. (1) mimics the availability of modern electronic communication between 
clusters and hence the possibility of some long-range inter-cluster 'coherence', while the fragmentation mimics the 
tendency for an individual cluster to scatter when suddenly sensing danger and hence lose its intra-cluster coherence 3 
(see Appendix A). The process of superconductivity enhancement and human conflict are therefore dynamically 
equivalent in that they can both be approximated by Eq. (1), and hence both produce identical exponentially 
truncated 2.5 power-law forms. We stress that we are not saying that these systems are physically identical, since 
obviously both the objects involved and the origins of their correlations are very different - but they do behave 
dynamically as 2/ they were both the same system (i.e. they both follow Eq. (1)). 

An immediate practical consequence is that learning how to 'protect' a population under attack from an insurgent 
or terrorist threat, and learning how to 'protect' a fragile SC state from decoherence and thus raise T c , become 
inter-related problems at this Master Equation level. Hence insights from one might help the other. This moves 
the Complexity debate forward from the issue of whether systems are physically similar or not, to whether they 
are dynamically equivalent or not at some given level of approximation (e.g. Master Equations for cluster sizes). 
To illustrate that it is not just SC systems that can belong to a given complexity universality class, we note that 
fragmentation in low-dimensional Bose disordered systems was recently identified as a novel mechanism for detecting 
a superfluid-insulator quantum phase transition^. 

This model (Eq. (1)) also offers an alternative explanation for a variety of other complex phenomena which have 
been found to exhibit a robust 2.5 power-law. Gabaix et a/P found a common power-law distribution for individual 
transaction sizes with a = 2.5 ±0.1, for the London Stock Exchange, the NYSE, and the Paris Bourse. Interpreting 
N as the average aggregate demand for stocks, this demand N gets shaped into a distribution of demand 'clusters' 
representing potential orders of a given size s. Since it is reasonable to expect orders to be realized at random, 
the distribution of individual transaction sizes is proportional to the distribution of clusters of potential orders - 
hence a = 2.5. Similarly, RichardsorP concluded that the distribution of approximately 10 3 gangs in Chicago, and 
in Manchoukuo in 1935, separately followed a truncated power-law with a ~ 2.3. Interpreting N as the number of 
potential gang members in each case, with each comprising a transient soup of clusters which tend to combine or 
fragment over time, yields a = 2.5. In a similar way, the robust time-dependence of a power-law with a « 2.4 in a 
recent New York garment industry stud}^can be reinterpreted as a repartitioning of trading interactions, with multi- 
component clusters continually being built up as part of common jobs (i.e. coalescence) and then dissolving upon 
completion (i.e. fragmentation). For collections of TV neurons^, we can imagine a dynamical coalescence- fragmentation 
grouping process in which groups of neurons become synchronized, and then this synchronization ultimately fragments. 
(Members of the same group need not be physically adjacent to each other). When an entire group fires, it creates a 
measurable activity equal to 5 the group size s. Hence the resulting activity distribution will follow a new power-law 
given by s x p(s). The resulting power-law exponent (a — 1) = 1.5, which is exactly the famous empirical 3/2 valued 
We note that although competing theories exist for many of these applications (i.e. Chialvo 5 for neural dynamics and 
Gabaix et al? for markets), we know of no other single mechanism which is simultaneously physically plausible for 
each application area and which can explain the mysterious finding of Fratini et a/.Hl. 

Additionally, the model yields several concrete predictions for superconductivity in the cuprates, and similar quan- 
tum phenomena, which we hope will stimulate future empirical investigation. First, it should be possible to manipulate 
host material properties such as to vary ^f rag and/or v coa \ and hence alter the cutoff value sq in Eq. (2). In particu- 
lar, if sq is measured across a wide range of samples, it should be possible to infer best-fit values for the dynamical 
quantities ^f rag and ^ coa i which directly affect the SC state's coherence, and which might have been impossible to 
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estimate by other means. Second, in the limit where the spatial extent of any potential phase coherence (i.e. cluster 
coalescence) is small, the model effectively becomes a low dimensional percolation problem in which the truncated 
power-law form now appears only at the percolation threshold pj^. As we move from the infinite dimensional (i.e. 
long-range interaction) model in Eq. (1) toward a two-dimensional limit with only nearest neighbor interactions, 
a —> 2 at pp**. This behavior has been analyzed by Laibowitz et alW* using a similar two-dimensional model, to 
account for the observed cluster statistics in thin films near p c . Away from p C: the distribution then completely loses 
its power-law characteristics, which is exactly what we have observed in the transition from modern insurgent wars 
(where long-range interactions are possible through communication devices and the Internet) to older wars which re- 
semble fights on a grid and in which interactions and information transfer are effectively limited to nearest neighbors. 
It would be fascinating to explore this same transition in La2Cu04 +2y (or an analogous SC system) by manipulating 
the fractal distribution of ordered i-Os during the annealing process. Third, it has been recently established that the 
model's temporal evolution is able to reproduce the large fluctuations associated with 28 human contagion in financial 
markets (currency trading on the one minute timescale), biological systems (cold transmission in schools), and social 
systems (YouTube downloads). It would be interesting to see if the SC fluctuations near T c are also well described 
by the stochastic numerical simulations of the model's temporal dynamics, and whether the averaged behavior of 
individual groups in the SC state near T c ends up following the analytic Master equations. Applying ultrafast optical 
pump-probe measurements on the SC state, it should be possible to track the fragmentation process and associated 
bursty behavior as a large phase-coherent cluster fragments and smaller clusters take over. Fourth, our finding raises 
the possibility that other types of quantum phase transition might also benefit from a space which better matches 
the possible fractal form of the underlying quantum many-body state. In particular, the survival of other exotic, yet 
fragile, quantum many-body states might be similarly enhanced by incorporating a protective 2.5 power-law 'box' 
within the host material microstructure. Such states may include highly nontrivial quantum correlations beyond 
two-body entanglement 29 . Indeed, a measurable shift in transition temperatures could be a novel, noninvasive way of 
probing such multipartite entanglement! 2 ^^. 



Appendix A 

Here we present the derivation of the exponentially cutoff 2.5 power-law. Analysis of a simpler version of Eq. (1), 
was completed by d'Hulst and Rodgers 23 , and real- world applications have focused on financial markets - however 
this is the first example of an application to superconductivity and derivation with general values ^f rag and v coa \. At 
each timestep, the internal coherence of a population of N objects (which we refer to as an 'agents' to acknowledge 
possible application to human systems) comprises a heterogenous soup of clusters. Within each cluster, the component 
objects have a strong intra-cluster coherence. Between clusters, the inter-cluster coherence is weak. An agent i is 
then picked at random - or equivalently, a cluster is randomly selected with probability proportional to size. Let S{ be 
the size of the cluster to which this agent belongs. With probability ^f rag , the coherence of a given cluster fragments 
completely into Si clusters of size one. If it doesn't fragment, a second cluster is randomly selected with probability 
again proportional to size - or equivalently, another agent j is picked at random. With probability ^ C oab the two 
groups then coalesce (or develop a common 'coherence' in the case of the superconducting application). The Master 
Equation is as follows: 

rj S— 1 ~ OO 

on^ _ I/coal y^, ■ / js Vfr a gsn s Zu coaX sn s ^ 

k=l k=l 

~W = ^£*V-^I> fc . (A2) 

k=2 * k=l 



Note here we make an approximation that TV — >• oo. The terms on the right hand side of Eq. (Al) represent all the 
ways in which n s can change. In the equilibrium state: 

1 - v s_1 

sn R = - - mg , „ krik(s — k)n s -k , s > 2 , (A3) 



(^frag + 2^ coa i)7V ^ 

^frag 
2(1 — Z/f ra g) 



m = ^~^n k . (A4) 



k=2 
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Consider 

CO CO 

G \.y\ = Yl knkyk = n iy + Yl knkyk = ni y + ' ( A5 ) 



A;=0 A: =2 



where y is a parameter and g[y] governs the cluster size distribution for k > 2. Multiplying Eq. (A3) by y s and 
then summing over s from 2 to oo, yields: 



i.e. 



5b] 2 - - 2n iy ) g[y] + njy 2 = . (A7) 

V ^coal / 



From Eq. (A5), g[l] = G[l] — n\. Substituting this into Eq. (A7) and setting y = 1, we can solve for g[l] 

g[1] = N (A8) 

^frag ~r ^^coal 



Hence 



m = N — g[l] = V{iae + i/coal AT . (A9) 

^frag H~ 2z/ CO al 



Substituting this into Eq. (A7) yields 



g[y]2 _ ( ^ ag + 2, coal jv _ 2iV(, frag + , coal ) \ + (iV(, frag + , coal ))^ 2 = Q (AiQ) 

\ ^coal ^frag + 2^ coa i / ( ^frag + 2 ^ coa i ) 

We can solve this quadratic for g[y) 



r i _ (^frag + 2^ coa i)7V / _ 4(^f rag + ^coal)^coal ^ _2 \\— ^ fra S + ^coal)^coal \ (All) 
4^ coa i I (^frag + 2^ coa i) 2 V (^frag + 2^ frag ) 2 i 



which can be easily expanded 



r i _ Qfrag + 2^ coa i)A/" ^ (2fc ~ 3)!! / 4Q frag + \* ( 

m ~ 2^ coal ^ (2fc)!! V (^fra g + 2^ coal ) 2 V " 1 j 



Comparing with the definition of in Eq. (|A5[) shows that 



We now employ Stirling's series 



^frag + 2^coal (gg 3)!! / 4Q frag + ^coaQ^coal V 

2^ coal s(2s)!! V Kag + 2^ coal ) 2 y 1 j 



ln[s!] = -ln[2n] + ( s + - ) ln[a] - s + — - ... . (A14) 



Hence for s > 2, we find 

'Ofrag + 2^ coa i)e 2 \ / 4(>f rag + ZWlHoal V ~ l) 2s ~ 3 / 2 



2 3 / 2 V^^coal A Kag + 2^ coal ) 2 / S 2 ^ 1 

which implies that 



N , (A15) 



^coal 

Qfrag + ^coal) 5 \ 5 /2 
Qfrag + 2^ coa i) 2 
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In the limit s ^> 1, this is formally equivalent to saying that 

n s ~ ex.p(-s/s )s~ 5/2 (A17) 

where 



so 



j n / 4(^f rag + ^coaQ^coaA 
(^frag + 2^coal) 2 / 



(A18) 



For large cluster sizes (i.e. large s such that s ~ O(N)) the power law behaviour is masked by the exponential 
function. The equilibrium state for the distribution of cluster sizes can therefore be considered a power-law with 
exponent a ~ 2.5, together with an exponential cut-off. The exponent a can be manipulated by implementing 
suitably chosen microscopic rules, as discussed by us and others elsewhere. 

In the human context, the fact that the interactions are effectively distance-independent as far as Eq. (1) is 
concerned, captures the fact that we wish to model systems where messages can be transmitted over arbitrary 
distances (e.g. modern human communications). Bird calls and chimpanzee interactions in complex tree canopy 
structures can also mimic this setup, as may the increasingly longer-range awareness that arises in larger animal, fish, 
bird and insect groups. In a human/biological context, a justification for choosing a cluster with a probability which 
is proportional to its size, is as follows: a cluster with more members has more chances of initiating an event. In 
the materials context, it will have a higher chance of feeling an interaction from a dephasing mechanism event (e.g 
impurity or phonon) and hence may fragment more readily. It will also be more likely to find members of another 
cluster more frequently, and hence be able to synchronize with them - thereby synchronizing the two clusters. It 
is well documented that clusters of living objects (e.g. animals, people) may suddenly scatter in all directions (i.e. 
complete fragmentation as in Eq. (1)) when its members sense danger, simply out of fear or in order to confuse a 
predator. Such fleeing behavior was discussed at length in the classic 1970 work 'Protean Defence by Prey Animals' 
by D. A. Humphries and P.M. Driver, Oecologia (Berl.) 5, 285-302 (1970). Clusters of inanimate objects such as 
doubly-ionized Argon atoms and animal Hox genes, also exhibit complete fragmentation. 
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